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As a result of a new improved fit to old bubble chamber data of the dominant axial C5 nucleon- 

to-Delta form factor, and due to the relevance of this form factor for neutrino induced coherent pion 

production, we re-evaluate our model predictions in Phys. Rev. D 79, 013002 (2009) for different 

observables of the latter reaction. Central values for the total cross sections increase by 20%~30%, 

^— V , while differential cross sections do not change their shape appreciably. Furthermore, we also compute 

fvj ■ the uncertainties on total, differential and flux averaged cross sections induced by the errors in the 

determination of C5 . Our new results turn out to be compatible within about la with the former 

^ ] ones. Finally, we stress the existing tension between the recent experimental determination of the 

^j ' CT(NCc'ohI-''') ^^^^^ by the SciBooNE Collaboration and the theoretical predictions. 
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I. INTRODUCTION 

Experimental analyses of neutrino induced coherent pion production generally rely on the Rein-Sehgal (RS) 
model [l|, [4) which is based on the partial conservation of the axial current (PCAC) hypothesis. In the RS model 
the pion-nucleus coherent cross section is written in terms of the pion-nucleon elastic cross section by means of ap- 
proximations that are valid for high neutrino energies and small values of the nucleus momentum transfer square and 
of the lepton momentum transfer square {q^). As pointed out in Refs. [3|, 0]) those approximations are less reliable 
for neutrino energies below/around 1 GeV, light nuclei, like carbon or oxygen, and finite values of q^ . These are the 
energies and targets used in present and forthcoming neutrino oscillation experiments. 

There are other approaches to coherent production that do not rely on PCAC but on microscopic models for pion 
production at the nucleon level [3|, I5l4l3| . The dominant contribution to the elementary amplitude at low energies is 
given by the A-pole mechanism (A excitation and its subsequent decay into iiN). Medium effects on the A mass and 
width, final pion distortion, as well as nonlocalities in the pion momentum, are very important and are taken into 
account in microscopic calculations. Similarly to PCAC models, the process is dominated by the axial part of the 
weak current and it is thus very sensitive to nucleon-to-Delta axial form factors. 

Very recently, the role of nonlocalities in the A momentum has also been investigated |13 - [l5| . In Ref. \\M, it is 
claimed that their neglect, the so called, local approximation, leads to an overestimate of the coherent production 
cross section that can be as large as a factor of 2 for neutrino energies of 500 MeV. Similar results were obtained in 
Ref. [15|. Final pion distortion and in medium modifications of the A properties were not considered, and it was not 
clear whether those approximations could not affect the results. Final pion distortion and in medium modification of 
the A properties were included in Refs. [l3;ll3l) were nonlocal effects on the A momentum were incorporated in the 
A self-energy in the first-order approximation. They also observe a large reduction in the total cross section due to 
the nonlocal aspects of the A propagation in the medium. However, as claimed by the authors of Ref. [Tjl, this that 
not mean that earlier microscopic calculations [3l.[5l-[ll| are wrong, as there A nonlocal effects are taken into account 
in an effective way through the in medium modification of the A properties which were fitted to observables. 

In the model we developed in Ref. Q, the A was treated in the local approximation. However, the modifications 
of the A in medium properties are such that similar models give a good reproduction of pionic atoms and tt— nucleus 



scattering [l6|, [iTf, pion photoproduction [l8|, pion electroproduction [l9|, (•^He,t) |20|] and elastic a— proton J21| 
reactions. We share the claim of the authors of Ref. [l^ and we believe this treatment of the A, where certainly 
non-local effects are being effectively (partially) taken into account, is also adequate for neutrino induced reactions. 
Nevertheless, this interesting issue deserves future investigations. 

Our model in Ref. Q is based on a microscopic model at the nucleon level, described in detail in Ref. [22;], that, 
besides the dominant A pole contribution, takes into account background terms required by chiral symmetry. As a 
result of the inclusion of background terms, we had to re adjust the strength of the dominant A pole contribution. 
The least known ingredients of the model are the axial nucleon-to-A transition form factors of which C^, not only 
gives the largest contribution, but it also controls all other axial form factors if one assumes Adler's model [231 that 

gives Ct{q^) = -^^1^, C^(g^) = 0, and PCAC is used to obtain Ci{q^) = C^{q^)-J^. This strongly suggested 
to us the readjustment of C5 to the experimental data. 

Information on pion production off the nucleon comes mainly from two bubble chamber experiments, ANL [24l426l | 

and BNL [23, |2a]. Assuming, as proposed in Ref. [231, the q^ dependence C^(g^) = .,_ 2/M'i -,2 (1 — o.Sa )~"^, we 

fitted in Ref. |22| the flux-averaged Vf^p -^ ii~p'K~^ ANL ^^-differential cross section for pion-nucleon invariant masses 
W < 1.4 GeVJll ^ obtaining C^(0) = 0.867 ± 0.075 and Maa = 0.985 ± 0.082, with a Gaussian correlation 
coefficient r ~ —0.85 and a y^ /dof = 0.4. The fitted axial mass was in good agreement with the estimates of about 
0.95 GeV and 0.84 GeV given in Refs. [H H^]. On the other hand, the C^ifi) value is some 30% smaller than the 
prediction obtained from the off diagonal Goldberger-Treiman relation (GTR) that gives C^(0)|g^^ = 1.2. C^(0) is 
not constrained by chiral perturbation theory (^PT) and lattice calculations are still not conclusive about the size of 
possible violations of the GTR. For instance, though values for C^(0) as low as 0.9 can be inferred in the chiral limit 
from the results of Ref. [3l[, they also predict C^(0)/ C's^CO)!^.^^ to be greater than one. 

Recently, two re-analysis have been carried out trying to make compatible the GTR prediction for C^{Q) and ANL 
data. In Ref. (S^], C'^(O) is kept to its GTR value and three additional parameters, that control the C^{q^) fall off, 
are fitted to the ANL data. Although ANL data are well reproduced, we find the outcome in [S^] to be unphysical, 
as it provides a quite pronounced g^— dependence giving rise to a too large axial transition radius of around 1.4 fm 
(further details are discussed in (3511 ^ ■ 

A second re-analysis [3^ brings in the discussion two interesting points. First that both ANL and BNL data were 
measured in deuterium, and second, the uncertainties in the neutrino fiux normalization. It is claimed in Ref. [33f, 
that the latter could be responsible for BNL total cross sections being systematically larger than ANL ones. In 



Ref. [33| the authors do a combined best fit to the ANL and BNL data including deuteron effects, which they evaluate 
as in Ref. [SJj, and flux normalization uncertainties, treated as systematic errors, and taken to be 20% for ANL data 
and 10% for BNL data. With a pure dipole dependence for C^, they found C^(0) = 1.19 ± 0.08, in agreement with 
the GTR estimate. 

The works in Refs. [321 133| consider only the A pole mechanism but ignore the sizable non-resonant contributions 
which are of special relevance for neutrino energies below 1 GeV. When background terms are considered, the tension 
between ANL data and the GTR prediction for C^(0) substantially increases as the results in Ref. [221 clearly shows. 

In our work in Ref. |^] we have performed a fit to both ANL and BNL data in which: i) We have included the 
BNL total Vf^p — > /i^p7r+ cross section measurements of Ref. [23]. We have just included the three lowest neutrino 
energies: 0.65, 0.9 and 1.1 GeV, since there is no cut in the outgoing pion-nucleon invariant mass in the BNL data, 
and we want to avoid heavier resonances from playing a significant role. We have not used the BNL measurement of 
the g^— differential cross section, since it lacked an absolute normalization, ii) We have taken into account deuteron 
effects, in) the uncertainties in the neutrino flux normalizations, 20% for ANL and 10% for BNL data, are treated 
as fully correlated systematic errors, improving thus the simpler treatment adopted in Ref. [33|, and finally iv) in 
some fits, we have relaxed Adler's model constraints, in order to extract some direct information on C^4^{0). For 

simphcity we took C^{q^) = ,^_ 2°;j^a -.2 ■ As in Ref. [33, the consideration of BNL data and flux uncertainties 

increased the value of C^(0) by about 9%, while strongly reduced the statistical correlations between C^(0) and 
Ma\- The inclusion of background terms reduced C^(0) by about 13%, and deuteron effects increased it by about 
5%, consistently with the results of [23 and [3J, Isj], respectively. Fitted data was quite insensitive to C^^{Q). 

In our most robust flt in Ref. ^g] we used Adler's constraints, and we obtained C^(0) — 1.00 ± 0.11, M^a = 
0.93 ± 0.07 GeV, with a small Gaussian correlation coefficient r — —0.06 and a x^ /dof — 0.32. This violation of the 
GTR is about 15%, and it is smaller than that suggested in Ref. [23, though it is definitely greater than that claimed 
in Ref. [3j|, mostly because in Ref. (331 background terms were not considered. However, the GTR value and the 
C^(0) above differ in less than 2(t, and the discrepancy is even smaller if Adler's constraints are removed. 

These new results are quite relevant for the neutrino induced coherent pion production in nuclei which is a low of 
dominated reaction. Background term contributions to coherent production largely cancel for symmetric nuclei [3| 
making the A pole mechanism the unique contribution. As the process is dominated by the axial part of the weak 
current, it is very sensitive to C^(0). Thus, we would expect the results in Ref. [3], based in the determination of 
C^(0) of Ref. [22|, to underestimate cross sections by some 30%. In this work we re-evaluate different pion coherent 
production observables using our model of Ref. [3| but with the new parameterization and results for C5 obtained in 
Ref. 36]. As the correlation coefficient is small in this case, we shall treat the theoretical errors that derive from the 
uncertainties in C^(0) and Maa as independent and we shall add them in quadratures. 



II. NEW RESULTS 

We start by showing in the left panel of Fig. [1] the differential cross section with respect to the £'7r(l — cos6'7r) 
variable for neutral current (NC) coherent 7r° production on carbon. E^^ is the pion energy in the laboratory frame 
(LAB) while Ot^ is the LAB angle between the pion and the incoming neutrino. The shapes of the distributions are 
completely similar to the ones we obtained in Ref. 'S], but the absolute values increase by some 20%^^ 30% depending 
on the neutrino energy. This is generally true for other differential cross sections that we do not show here. For the 
distribution convoluted with the MiniBooNE flux we flnd an increase in the total cross section of about 29%. 

In Table HI we show our new predictions for, both NC and charged current (CC) processes, for the K2K [37] and 
MiniBooNE [331 ^^^^ averaged cross sections as well as for the future T2K experiment. In the middle and right panels 
of Fig. [Tl we show some results for T2K and MiniBooNE experiments. In all cases, the flux is normalized to one. As 
in Ref. [3|, and since we neglect all resonances above the A(1232), we have set up a maximum neutrino energy (i?^a,x) 

in the flux convolution, approximating the convoluted cross section by tr « J^"^""" dE(f>'^{E)(T{E)/ J^,""''' dE(j)'^{E), 

where we fixed the upper limit in the integration to -Emax — 1.45 GeV and 1.34 GeV for CC and NC v^/i^^ driven 
processes, respectively. E^^,„ is the lower flux limit. For the K2K case a threshold of 450 MeV for muon momentum 
is also implemented [33l and we can to go up to i?^£;.^^^=1.8 GeV. We cover about 90% of the total flux in most of 
the cases. For the T2K antineutrino flux, we cover just about 65%, and therefore our results are less reliable. 

Our central value cross sections increase by some 23%~30%, while the errors associated to the uncertainties in the 
C^(0) and Maa determination are of the order of 21%. Our new results are thus compatible with former ones in 
Ref. [3 within la. Our prediction for the K2K experiment lies more than la below the K2K upper bound, while 
we still predict an NC MiniBooNE cross section notably smaller than that given in the PhD thesis of J.L. Raaf J38|. 
Note however, that the MiniBooNE Collaboration has not given an official value for the total coherent cross section 
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TABLE I. NC/CC v^ and Vf^ coherent pion production total cross sections, with errors, for K2K, MiniBooNE and T2K experiments. 
In the case of CC K2K, the experimental threshold for the muon momentum |fc,i| > 450 MeV is taken into account. Details on the flux 
convolution are compiled in the last three columns. 



yet, and only the ratio coherent/ (coherent+incoherent) has been presented [30|. For the future T2K experiment, we 
now get cross sections of the order 2.4— 3.2xl0~''°cni^ in carbon and about 2.9—3.8 xlO~''*'cm^ in oxygen. 

In Fig.[2]we show new v^/v CC and NC coherent pion production total cross sections off carbon and oxygen targets. 
As in Ref. [3|, we observe sizable corrections to the approximate relation ace , ~ 2(Tnc for these two isoscalar nuclei 
in the whole range of v/v energies examined. As pointed out in Refs. [2|, |40|, this is greatly due to the finite muon 
mass, and thus the deviations are dramatic at low neutrino energies. In any case, these corrections can not account 
for the apparent incompatibility among the CC K2K cross section and the NC value quoted in Ref. |38| . 

The SciBooNE Collaboration has just reported a measurement of NC ■n'^ production on carbon bj^ a v^j, beam with 
average energy 0.8 GeV [4l[ 
that iltcccoh^ 



T(NCcoli7rO) 



SciBooNE 



Based on previous measurements of CC coherent 7r+ production [42|, they conclude 
= 0.14^Q2g. This result can not be accommodated within our model, or any other present 



theoretical model, neither microscopic pj|-[Tll. Il3j nor PCAC based [ilHjSj,!^- Theoretically, this ratio cannot be 
much smaller than 1.4-1.6. For instance, for a carbon target and for a neutrino energy of 0.8 GeV we find a value of 
1.45 ± 0.03 for that ratio, ten times bigger that the value given by the SciBooNE Collaboration. From the v^ +^^ C 
CC and NC MiniBooNE convoluted results shown in TablelH we obtain 1.46 ± 0.03. We believe this huge discrepancy 
with the SciBooNE result stems form the use in Ref. 4Jj] of the RS model to estimate the ratio between NC coherent 
7r° production and the total CC pion production. As clearly shown in Refs. [3, 43, the RS model is not appropriate to 
describe coherent pion production in the low energy regime of interest for the SciBooNE experiment. 
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FIG. 1. Left panel: Laboratory E-^(l — cos9,r) distribution for the v C — > u Ctt reaction, at MiniBooNE energies. We also show, 
with the corresponding error band, the distribution convoluted with the u^ MiniBooNE flux. We display the MiniBooNE published 
histogram, taken from the right panel of Fig. 3 in Ref. [3S(I . conveniently scaled down so that their first bin matches our result at 
E„{1 — cos 9^) = 0.005 GeV. Middle and right panels: CC and NC coherent pion production cross sections in carbon (dots). We also 
show (solid lines) predictions multiplied by the T2K (middle) and MiniBooNE (right) Ufj, energy spectra. The dashed curves curves stand 
for the T2K and MiniBooNE u^ fluxes normalized to one. Error bands are shown for our results. 
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FIG. 2. Ufi/Ufi CC and u/u NC coherent pion production from a carbon target (left/middle panel) and i^^ CC and u NC coherent pion 
production from an oxygen target (right panel) as a function of the neutrino/antincutrino energy. Error bands are shown. 
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